Optimal non-linear passage through a quantum critical point 
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We analyze the problem of optimal adiabatic passage through a quantum critical point. We show 
that to minimize the number of defects the tuning parameter should be changed as a power-law 
in time. The optimal power is proportional to the logarithm of the total passage time multiplied 
by universal critical exponents characterizing the phase transition. We support our results by the 
general scaling analysis and by explicit calculations for the transverse field Ising model. 



The adiabatic quantum computation is a very appeal- 
ing concept for simulating certain class of problems [l|, [3] . 
The main idea is that the system is initially prepared in 
a ground state of some simple Hamiltonian and then the 
Hamiltonian is adiabatically changed to some nontriv- 
ial regime. If one changes the Hamiltonian sufficiently 
slowly then the system remains in the ground state and 
this new ground state may encode nontrivial information. 
For example the adiabatic evolution can be a necessary 
condition for realizing topological quantum computation 
in cold atom systems such as p-wave superfluids ^s', '3| . In 
order to get to the nontrivial topological state one needs 
to cross a second order phase transition between gapp- 
less and gapped px + ipy phases [1] . Indeed it is generally 
much easier to cool the system in the regime where topo- 
logical excitations do not exist and then transfer the sys- 
tem to the nontrivial phase. Similar considerations can 
be applied to proposals based on the Kitaev model 5J. 
We note that this idea can be used in a broader sense than 
quantum computing alone. For example, one can use the 
same algorithm to study ground states of strongly corre- 
lated systems. The adiabatic evolution is routinely used 
in cold atom systems to load them into optical lattices, 
change interaction strength, etc. Q. 

Of course ideally it is good to prepare the system ini- 
tially in the ground state and keep the evolution slow 
enough so the system is never excited However this 
can be a very stringent requirement. If a sufficiently large 
system is cooled in the gapless regime it is virtually im- 
possible to avoid any excitations. In addition the require- 
ment that the evolution is completely adiabatic might 
lead to extremely long computational times and make 
the whole idea of adiabatic quantum computation unre- 
alistic. At the same time the quantum computation can 
be robust to certain appropriate level of excitations and 
thus one should always look for a compromise between 
shortening the evolution time and tolerance to the num- 
ber of the excited states. Likewise if one is interested in 
understanding many-body states using adiabatic evolu- 
tion of the Hamiltonian one should compromise between 
the evolution time and the heating generated in the sys- 
tem due to non-adiabaticity. Usually the effects of non- 
adiabaticity are particularly strong in low dimensional 
systems 0] especially near singularities such as critical 
points. By now there is an extensive literature focusing 



on slow dynamics near critical points 
and even critical lines like in the Kitaev model [14!]. Typ- 
ically as a measure of non-adiabaticity one defines the 
density of produced defects during the quench. Such 
measure is appropriate for integrable models with well 
defined excitations. More generally one can look into the 
excess energy density added to the system during the 
quench and come to similar conclusions 

Gencrically if one crosses a phase transition, the tuning 
parameter g will change linearly in time near the critical 
point: g t/T, where T is the characteristic passage 
time. For convenience we assume that g vanishes right 
at the critical point. If one is interested in minimizing the 
heating in the system or the number of produced defects 
in the system then such a linear change might not be 
optimal. Indeed it is intuitively clear that one should slow 
down precisely near the critical point and perhaps change 
the tuning parameter as a power law in time: g ^ {t/Ty . 
Now one can ask the question: what is the optimal value 
of r? Is it determined by some universal properties of the 
phase transition? Or perhaps there is an optimal value 
of the power, e.g. r oo, which is good irrespective 
of the details of the critical point. The answer to these 
questions is the main purpose of the present letter. 

Let us first argue that there should be an optimal value 
of r. Indeed if r — > then the tuning parameter abruptly 
changes near the critical point, and we expect many ex- 
citations to be produced in the system. On the other 
hand if r — > oo then the tuning parameter changes very 
smoothly near the critical point but then changes very 
abruptly inducing a very strong heating in the system. 
If we want to find an optimal path at a given constraint 
that the total passage time T is fixed we clearly have 
to find a compromise between being as slow as possible 
right at the critical point and being not very fast when 
the tuning parameter is finite. 

Let us now be more specific and assume that 



g(i)«go|t/rrsign(t). 



(1) 



where go is the characteristic value of the tuning param- 
eter in the final state. Note that the time t can change 
either within [— T, T] corresponding to the crossing of 
the critical point or within [0, T] corresponding to start- 
ing right at the critical point and changing the tuning 
parameter from zero to the final value go. The latter 
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regime is realized, for example, if one loads interacting 
one-dimensional Bose gas into an optical lattice A 
generic quantum phase transition at 5 = is character- 
ized by critical exponents 17[. The exponents, which are 
relevant to us are v and z. The correlation length ex- 
ponent V characterizes the divergence of the correlation 
length ^ at small g: f ~ ^ while the dynamical critical 
exponent z defines the relation between a characteristic 
spectral energy A and the correlation length: A ^ 
(see Ref. 17| for more details). Instead of T it is con- 
venient to introduce a dimensionless parameter of adia- 
baticity, J = l/(TAo), where Aq corresponds to g = go- 
For gapped systems Aq is related to the lowest excitation 
energy in the system, while for gapless systems Aq is the 
energy scale at which the spectrum qualitatively changes 
from the low energy to the high energy behavior 17[ • Our 
main conclusion is that in the adiabatic limit 5 <IC 1 the 
optimal power Topt depends only on (5 and the universal 
exponents characterizing the phase transition: 
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ln(C/<5) 



(2) 



where C is a non-universal constant of the order of unity, 
which depends on the details of the phase transition and 
precise choice of Aq in the definition of the parameter b. 

The result ([2]) is quite remarkable. It implies that the 
optimal power (modestly logarithmically) depends on the 
total passage time. It approaches infinity only for in- 
finitesimally slow processes. If the passage time is large 
but finite then ropt is usually a relatively small power of 
the order of unity. Interestingly dimensionality does not 
enter the result ([2]). Similarly, we get the result for the 
density of the produced defects if one uses the optimal 
adiabatic passage: 
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ln(C/,5) 



(3) 



Here d is the dimensionality of the system. In dimensions 
higher than one these results immediately generalize to 
the situation when instead of critical point one has a 
critical surface of dimension d\ 



Kitaev model. Then in Eq 
d — dii instead of d 
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E.g. d|| = 1 for the 
([3]) one has to use d\^ = 
The scaling of the density of the 
produced defects for optimal passage is quite different 
from what one would get for a generic linear in time 
crossing of the critical point [1]: ^ §dv/(zij+i) ^ j^- 
is easy to check that: 
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< 1. 



(4) 



so that optimal passage gives a significant improvement 
over the linear one. We note that the constant C can 
always be reabsorbed in the precise definition of 5. To 



simplify the notations we will assume this is the case and 
we skip C in the expressions below. 

We note that the optimal path docs not have to be a 
power law function at t — > as it is assumed in Eq.(IT]). 
However, the power law behavior is very universal and 
is expected to be valid close to the critical point. The 
specific details of the optimal functional dependence on t 
away from the critical point are likely to be non-universal. 
However, since the defects are mainly produced at g{t) 
close to zero, we do not expect that our results will be 
qualitatively affected by the deviations from the power 
law. In the rest of the paper we will sketch the derivation 
of our main result and consider a specific example of the 
transverse field Ising model. 

To calculate the number of the defects produced dur- 
ing the process we will use the Fermi Golden rule anal- 
ysis in the adiabaticity parameter 6. As it was shown 
in Refs. 0, Q in the case of linear quench unless the 
non-adiabtic regime is realized 01 this analysis gives the 
correct dependence of the density of the defects on 5 but 
does not reproduce the prefactor correctly. For exam- 
ple, the Fermi Golden rule calculation overestimates the 
prefactor by approximately two times for the transverse 
Ising model. Technically the mistake comes from the fact 
that at very small energies the process is nonadiabatic 
and the Fermi Golden rule is not precise. However, one 
can show, that unless the system is in the non-adiabatic 
regime the mistake reduces to a mere factor, which is al- 
ways bounded between one and another constant of the 
order of one. This result is true for any choice of time 
dependence of the tuning parameter g thus the general 
scaling results based on the Fermi Golden rule analy- 
sis are correct. We will illustrate the difference between 
Fermi Golden rule and exact approaches below using a 
particular example of the transverse field Ising model. 
As we will show there is only a minor difference between 
them. 

Within the Fermi Golden rule the density of the excited 
states is given by the following expression: Q 



dt{p,t\dt\0,t)efo''''i^-^''^-'^°^''^^ 



(5) 



where the last summation is over the d-dimensional phase 
space. Here, \p, t) is the instantaneous eigen-state of the 
Hamiltonian describing the system, and Ljp{t) is the exci- 
tation energy. We can next use general scaling arguments 
(see also Ref. [8] ) to find the dependence of n on the rate 
5 and the exponent r. Thus for the excitation spectrum 
we can use 

- = f'fiNp') = f-f mr^'/pi , (6) 

where f{x) is a scaling function satisfying f{x) ~ 1 at 
X <C 1. For gapped systems we should also have f{x) ^ x 
at a; 3> 1. The latter condition implies that Up — loq ^ A 
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for p — >■ 0. If the system is gapless then the asymptotics 
of / will be somewhat different and our analysis should 
be modified accordingly. 

Dependence of the gap on time, A = ji^j''^'' suggests 
rescaling the variables 



T=- , q^p[S{rzi^ + l)] — +1 , (7) 

which leads to the following expression for the phase de- 
fined as (fiq — /Q((jjp(i') — uJo{t'))dt': 

V,{x) = f 1 + — ) <?^(^+^> r dyy^-' f{y), (8) 
\ rzv/ Jq 

where x{t) — sign(T)|T|''^''/q^. Similarly one can rescale 
the matrix element in Eq.Q, dt{p, t\dt\0, t) — dxG{x), in 
which G{x) « const at |a;| <C 1. 

The main contribution to the number of defects in 
Eq. ([5]) can be extracted from the behavior of the phase 
factor in Eq.®, which is controlled by the asymptotics 
of the scaling function f{x). At \x\ ^ 1, the phase can 
be approximated by </?g(a;) ~ q^'^^^'^'> [rzv + l)x^/'^'"^ ^ 
while at 3> 1 it is (pq{x) « (g^a;)^^^^) approach- 
ing a linear function at r ^ 1. The latter asymptotics 
nearly independent of the loading power corresponds to 
the dominant contribution of the excited states to the 
number of defects. The characteristic rescaled momen- 
tum g^,, ~ 1 at which the transitions become suppressed 
provides an estimate for the width of excited energy 
states = (r5)'-'^/('-^'^+i). 

Combining the scalings above we find 

n ^ ^ (9) 

We note that the scaling of n with 5 for an arbitrary 
power r was recently obtained in Ref. [isj using similar 
methods. He we generalized that result to include the 
dependence of the scaling on power r, which is crucial for 
finding the optimum. Indeed it is this dependence, which 
encodes non-adiabatic transitions away from the critical 
point. Let us point that Eq. ([9]) agrees with our general 
expectations outlined in the beginning of the text. Thus 
at small r the density of excitations is large because the 
exponent of S is small. On the other hand for rzv 3> 1 the 
exponent of S saturates at a constant value and does not 
depend on r any more while the prefactor grows with r. 
Clearly for any given 6 there should be an optimal value 
of r. And indeed if one extremizes Eq. ^ with respect 
to r one finds that there is a minimum of n occurring at 

ToptS « e-'--'-^ (10) 

At (5 <C 1 this equation results in Eqs. ^ and 

To illustrate our results we will use a specific one- 
dimensional transverse field Ising model defined by the 
Hamiltonian: 

i i 
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FIG. 1: The prefactor A{r) in the number of defects Eq. (fTS)) 
as a function of the loading power for the passage through the 
critical point in the transversed field Ising model: the exact 
calculation (red solid line) and perturbative analysis (black 
solid line) of the adiabatic passage. Inset: prefactor A{r) for 
the loading. 

where af and af denote Pauli matrices. This model 
describes a second order phase transition between two 
gapped phases at g = [l7| (we assume that g < I). The 
critical indices corresponding to this phase transition are 
V = z = 1. We will consider two generic possibilities: "a 
passage", where one starts deeply in the one phase and 
changes g through the critical point and "loading" where 
one starts exactly at the quantum critical point g — 
and increases the tuning parameter to drive the system 
into one of the phases. The second possibility, for exam- 
ple, naturally arises if one loads commensurate bosons to 
a one-dimensional periodic potential [l^. The problem 
maps exactly to the transverse field Ising model with the 
amplitude of the lattice potential being proportional to 
1 - .9 ■ 

Using the Jordan- Wigner transformation one can map 
the Hamiltonian (jlip to the Hamiltonian of noninteract- 
ing fermions (l7j : 

^ = E (^I^'+i + 4+i4 + (9- l)clc. + h.c.) . (12) 

i 

This Hamiltonian can be diagonalized by standard Bo- 
goliubov transformation. Then the excitations corre- 
spond to the pairs of fermions moving with opposite mo- 
menta. Near the critical point the excitation energy of a 
pair and the matrix elements are given by 

LOp « 4Vp2+52, (13) 

and 

(P.5|9,|0,5)«^-^. (14) 
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The factor of four in the equation for the spectrum comes 
from the fact that the quasiparticles are excited in pairs 
with opposite momenta. The matrix element (|14p corre- 
sponds to the pair of excitations with momenta p and —p. 
We assume that g = go{t/TYs\gn[t) and t £ [-T, T]. We 
note that for large T the number of excitations does not 
depend on the upper limit so we will choose 50 = 1 
and then 5 = l/(rgo) = 1/^. Using Eqs. and 
([5|) it is straightforward to show that within the Fermi 
Golden Rule the expression for the density of excitations 
at small values of 5 assumes the form: 



[A{r)5 



r/(r+l) 



(15) 



where within the perturbation theory 



A{r) = {r + l) 



dq 
2^ 



dx 



1 + x- 



■ cos (Pq{x) 



(16) 



The phase (pqr{x) at r ^ 1 in Eq. (|16p is a function of a 
rescaled momentum q which reaches the order of unity at 
a characteristic momentum q^ = r' /'^''^^-'. By evaluating 
the integral in Eq. ([5]) one obtains the scaling with 
with the linear function Aper{r) ~ 0.039(r — 0.2) for the 
passage and Aper{r) « 0.02(r — 0.12) for the loading [3]. 
For the transverse field Ising model the density of the 
produced defects can be also found exactly since for each 
momentum one needs to solve an independent Landau- 
Zener problem [10]. The result of this calculation repro- 
duces the dependence ()15p with the following asymptoti- 
cal behavior of Aex{r) at large r: Aex{r) « 0.025(r — 1.3) 
for the passage and Aex{r) ~ 0.017(r — 1.3) for the load- 
ing. The linear scaling of A(r) with r is in complete 
agreement with general argument presented above if one 
sets z = 1/ = 1. 

The exact prefactors A{r) are compared in Fig[l]with 
approximate ones obtained from the pcrturbative calcu- 
lation in Eq. pB)) . In all cases A{r) linearly increases with 
r at large r. As we mentioned earlier the ratio of the ex- 
act and golden rule result is a number which only slightly 
depends on r and saturates at a constant independent of 
r at large r. 

Minimizing the density of excitations with respect to 
the power r, one arrives at the optimal power of loading: 
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IniC/S) 



(17) 



where for the passage constant C = 14.7 is extracted 
from the linear fit of the prefactor A^x {r) , and similarly 
C = 21.6 for the loading. It is a weak function of the 
dimensionless loading parameter S. This equation agrees 
with our general result 12]). The analytic expression for 
the optimal loading power gives the right asymptotic be- 
havior at small values of the adiabaticity parameter as 
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FIG. 2: The optimal number of excitation as a function of 
the adiabaticity parameter for the passage through the crit- 
ical point in the transversed field Ising model. The result of 
the exact numerical calculation (red solid line) is compared 
with analytic estimate (dashed blue line) and linear passage 
(black solid line). Inset: the optimal loading power as a func- 
tion of adiabaticity parameter for the passage through the 
critical point in the transversed field Ising model: the exact 
calculation (solid red line) and analytic estimate (dashed blue 
line) in Ea. ([T7|l . 



shown in Figj^l The dependence of the optimal number 
of excitations on the adiabaticity parameter is shown in 
Figl21 At small values of S, Eq.© is a good approxi- 
mation of the exact result. The optimal passage gives 
significantly less defects than the linear one. 

In conclusion we showed that there is a universal op- 
timum adiabatic path which minimizes the number of 
excitation produced in the system if the latter is slowly 
driven through a quantum phase transitions. This path 
is characterized by a power law dependence of the tun- 
ing parameter near the critical point, where the power 
is proportional to the logarithm of the total passage 
time multiplied by universal critical exponents charac- 
terizing the phase transition. Our results can be used 
both for optimization of adiabatic quantum computation 
algorithms and for minimizing heating in cold atom sys- 
tems in strongly correlated regimes. 
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